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ABSTRACT 

Israelian et al. (1999) showed that the stehar companion of the black-hole binary Nova Sco is polluted 
with material ejected in the supernova that accompanied the formation of the black- hole primary. Here 
we systematically investigate the implications of these observations for the black- hole formation process. 
Using a variety of supernova models, including both standard as well as hyper nova models (for different 
helium- star masses, explosion energies, and explosion geometries) and a simple model for the evolution of 
the binary and the pollution of the secondary, we show that most of the observed abundance anomalies 
can be explained for a large range of model parameters (apart from the abundance of Ti). The best 
models are obtained for He star masses of 10 to 16 M©, where spherical hypernova models are generally 
favoured over standard supernova ones. Aspherical hypernova models also produce acceptable fits, 
provided there is extensive lateral mixing. All models require substantial fallback and that the fallback 
material either reached the orbit of the secondary or was mixed efficiently with material that escaped. 
The black hole therefore formed in a two-step process, where the initial mass of the collapsed remnant 
was increased substantially by matter that fell back after the initial collapse. This may help to explain 
the high observed space velocity of Nova Sco either because of a neutrino-induced kick (if a neutron star 
was formed first) or by asymmetric mass ejection in an asymmetric supernova explosion. 

Subject headings: black holes — stars: binaries — stars: evolution — stars: individual (Nova Sco) — 
super novae: general 



1. INTRODUCTION 

X-ray Nova Sco 1994 (GRO J1655-40; hereafter Nova Sco) 
is one of the best-studied black-hole transient of recent 
years (e.g., Bailyn et al. 1995; Harmon et al. 1995; Hjellm- 
ing & Rupen 1995; Tingay et al. 1995; Orosz & Bailyn 
1997; van der Hooft et al. 1998; Shahbaz et al. 1999). It 
is a low-mass black-hole binary with an orbital period of 
2.61 d and relatively well-determined component masses. 
The most recent and most self-consistent analysis of the 
ellipsoidal light curves of the system by Beer & Podsi- 
adlowski (2001) has yielded masses of 5.4 ± 0.3 and 
1.45 ± 0.35 M0 for the black hole and the secondary, re- 
spectively, which we adopt in this study^. 
Nova Sco stands out among the black-hole transients be- 

^ The studies by Orosz Bailyn (1997) and van der Hooft et al. (1998) obtained somewhat higher masses for both components. However, 
both of these studies used a color excess, E{B — y) ~ 1.3, that is too large to be consistent with the observed colors of the secondary in Nova 
Sco, assuming that its spectral type is in the range of F2-F7 HI/IV (see Beer & Podsiadlowski 2001 for details). 



cause of an unusually high space velocity. Using the 7- 
velocity of Shahbaz et al. (1999) and the corrections for 
Galactic rotation of Brandt, Podsiadlowski & Sigurdsson 
(1995), Nova Sco's space velocity must exceed 106kms~^, 
a factor of a few larger than in any other low-mass black- 
hole transient. Brandt et al. (1995) concluded that the 
most likely explanation for the high space velocity is that 
the black hole formed in a two-stage process where the 
initial collapse led to the formation of a neutron star ac- 
companied by a substantial kick (Lyne & Lorimer 1994). 
The neutron star was subsequently converted into a black 
hole by accretion of matter that was not ejected in the 
supernova or a phase transition in the cooling compact 
object (Brown & Bethe 1994). (For a different view, see 
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Nelemans, Tauris & van den Heuvel [1999].) 
That the black hole formed in a supernova event was con- 
firmed by Israelian et al. (1999). Using high-resolution 
Echelle spectroscopy with the Keck Telescope, they 
showed that the atmosphere of the secondary was enriched 
by a factor of 6-10 in several o^-process elements (O, Mg, 
Si, S, Ti; see Table 1). Since some of these elements are al- 
most exclusively synthesized during a supernova explosion 
and cannot have been produced in a low-mass secondary, 
the secondary must somehow have been exposed to super- 
nova material that was ejected when the compact object 
in Nova Sco was formed. It should be noted that it is by 
no means required that the formation of a black hole is ac- 
companied by a supernova-like event. If a black hole forms 
promptly, i.e., on the dynamical timescale of the collaps- 
ing core, very few neutrinos can escape from the collapsing 
object (e.g., Gourgoulhon & Haensel 1993). Since delayed 
neutrino heating may be an essential feature in producing 
a successful ejection of the stellar envelope, i.e., a super- 
nova, a lack of neutrino emission may lead to a "failed" su- 
pernova in which the whole star collapses into a black hole. 
In this case, no pollution of the secondary is expected, nor 
is a supernova kick. The observations of Israelian et al. 
(1999) immediately rule out a simple prompt black- hole 
formation scenario for Nova Sco without an accompany- 
ing supernova and suggest that the black hole formed in a 
two-step process, where a neutron star may have formed 
first and was subsequently converted into a black hole by 
accretion. Alternatively, a black hole may have formed 
promptly, but subsequent accretion from a surrounding 
disk could have driven a jet-like explosion, as in the col- 
lapsar models of MacFadyen & Woosley (1999). In either 
case, the compact star could have received a substantial 
kick, in the first case due to an asymmetry in the neutrino 
emission (as may be the cause of neutron-star kicks), in 
the latter perhaps because of an asymmetry in the jets. 
The purpose of the present study is to explore the impli- 
cations of the observations of Israelian et al. (1999) for the 
formation of the black hole in Nova Sco in some detail. 
While Si is enriched by a factor of 8, Fe - surprisingly - 
is not. Since both of these elements are produced close to 
the mass cut above which matter is ejected in a success- 
ful supernova, these observations directly probe the region 
that is most crucial in determining whether a supernova 
is successful or fails. Indeed, the observations pose an 
immediate problem, since, in a standard supernova, the 
mass cut (Mcut — 1.5-2 M©; Thielemann, Nomoto, & 
Hashimoto 1996) is much smaller than the present mass 
of the black hole in Nova Sco. Figure la shows the com- 
position of the ejecta for a 16 Mq helium star for a stan- 
dard supernova (i.e., with a canonical explosion energy 
E'k = 1 X 10^^ erg; from Nakamura et al. 2001). Elements 
like S and Si, both enhanced significantly in the compan- 
ion of Nova Sco, are produced between the innermost 2.5 
and 3.5 Mq of the helium core, well below the final mass 
of the black hole. There are four possible ways by which 
some matter synthesized in this region can reach the sec- 
ondary and by which the mass of the black hole can be 
increased to the present value: (1) by fallback of material 



in the supernova explosion (Woosley & Weaver 1995), (2) 
by post-supernova mass transfer from the secondary, (3) 
as a result of mixing during the collapse phase and (4) as 
a result of a more energetic supernova (a hypernova with 
a larger mass cut). The composition of a hypernova model 
(with = 3 X 10^^ erg) is shown in Figure lb. It illus- 
trates how, for a more energetic hypernova, elements such 
as S and Si are produced much further out in the core of 
the helium star. Indeed it was this model that motivated 
Israelian et al. (1999) to first suggest a hypernova model 
for the formation of the black hole in Nova Sco. This and 
a possible connection to gamma-ray bursts was further de- 
veloped by Brown et al. (2000). These various solutions 
are, of course, not mutually exclusive, and we shall con- 
sider all of them in the following sections (see Fig. 2 for a 
schematic picture of the various cases considered). 

2. MODELING THE POLLUTION IN THE SECONDARY 

In order to examine whether the observed pollution of the 
secondary is consistent with the predictions of stellar and 
supernova nucleosynthesis, we have to follow the evolution 
of the binary through various evolutionary phases before 
and after the supernova and need to model the pollution 
of the secondary. In this section we present a very simple 
model that includes the main physical effects. Its purpose 
is to serve as a reference model with which we can dis- 
cuss the physical implications of the observed abundance 
anomalies. In § 3 we critically assess some of the assump- 
tions in this model and successively add physical realism 
to the model by first including mixing in the ejecta (§ 3.3, 
3.4) and then considering an aspherical explosion (§ 3.5). 
As the model becomes more realistic, the modeling uncer- 
tainties also increase. 

2.1. Description of the model without mixing (Case A) 

The calculation of the pollution of the secondary requires 
stellar models which have been subjected to explosive nu- 
clear burning in a supernova event. We take these from 
the library of models calculated by Nomoto et al. (1997, 
2001a, b). At the time of the explosion, the immediate su- 
pernova progenitor had to be a helium star (or Wolf-Rayet 
star) in order to fit into the tight binary orbit implied by 
the present orbital period of Nova Sco. We consider helium 
stars with initial helium-star masses, Mj^^, of 6, 8, 10 and 
16 M0. These correspond to main-sequence masses of ^ 
20, 25, 30, 40 M©, respectively. We consider two classes of 
explosion models: supernova models that have a standard 
supernova explosion energy of 1 foe (1 foe =10^^ erg) and 
hypernova models with explosion energies of 8 and 30 foe 
(for the 10 M0 model) and 30 foe (for the 16 model). 
The 6 Mq model is a model calculated for SN 1987A, while 
the 30 foe, 16 Mq model is appropriate for the prototype 
hypernova SN 1998bw (Iwamoto et al. 1998; Nakamura et 
al. 2001)2. 

Helium/ Wolf-Rayet stars are known to lose a substantial 
fraction of their envelopes in a stellar wind. To take this 
into account, we assume that the helium star has lost an 
amount AMh© = g (^He ~ ^bh) before the explosion, 

^ The term 'hypernova' was coined by Paczyhski (1998) as a model for gamma-ray bursters, linking them to the cataclysmic deaths of massive 
stars and the formation of black holes, as in the failed- supernova model of Woosley (1993) (for more details see MacFadyen & Woosley 1999). 
For the purposes of this study, we just define them as very energetic supernovae with energies > 10^^ erg, as has been used by Iwamoto et al. 
(1998) and Nomoto et al. (2001a,b). 
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where is the initial mass of the compact remnant 

(neutron star or black hole). We use the parameter g to 
vary the total amount of wind mass loss before the super- 
nova. 

At the time of the explosion, the masses of the primary and 
secondary are Mhg and , respectively. When the pri- 
mary collapses, it first forms a compact remnant of mass 
Mgfj. The rest of the envelope is assumed to be ejected 
initially, but part of it (Mfaiiback) will subsequently fall 
back, either because it did not achieve escape velocity or 
was pushed back by a reverse shock in the envelope (see 
Woosley & Weaver 1995). The fallback matter increases 
the mass of the compact remnant to (indeed, it may 

be this fallback that leads to the conversion of the com- 
pact remnant into a black hole). Figure 3a schematically 
illustrates the definition of these various mass parameters. 
In the simple model we consider first, we assume that 
all of the matter that falls back has moved beyond the 
position of the secondary (in § 3.2 we shall critically as- 
sess this assumption). Therefore, the secondary can be 
polluted twice with supernova material, first by all the 
material that is ejected and then by material that falls 
back. We assume that the fraction of matter that is cap- 
tured is given by the geometric fraction of the secondary 
([i?o/2ao]^ 0^ 0.01-0.03, where = O.SRq {M^ M^f-^ 
is the radius of the secondary and ao the initial orbital 
separation) times some efficiency factor /, where we as- 
sume different efficiency factors for matter that passes the 
secondary in the initial ejection (/ejection) and for matter 
that falls back (/fallback)- The efficiency factors can be 
much smaller than 1, for example, if the supernova leads 
to stripping of matter from the secondary (Marietta, Bur- 
rows, & Fryxell 2000 and § 3.1), or larger than 1 if grav- 
itational focusing is important. The latter requires that 
the relative velocity of the material is less than the escape 
velocity of the secondary and can plausibly only occur for 
fallback material. We also take into account the pollution 
of the secondary that has occurred before the supernova 
because of the capture of wind material by the secondary 
(where we assume a capture efficiency of 1). 
The matter that is captured by the secondary has a much 
larger mean molecular weight than the composition of the 
secondary, a relatively unevolved star at this stage. This 
is secularly unstable and leads to thermohaline mixing in 
the secondary (e.g., Kippenhahn, Ruschenplatt & Thomas 
1980). Since the time scale for thermohaline mixing is 
short compared to the evolutionary time scale of the sec- 
ondary, we assume that the material captured by the sec- 
ondary is completely mixed with the rest of the star after 
the supernova. 

In order to be able to follow the post-supernova evolution, 
we assume that the pre-supernova system is circular and 
that the supernova explosion is spherically symmetric in 
the frame of the primary. It is then straightforward to es- 
timate the post-supernova parameters of the system (we 
follow Brandt & Podsiadlowski 1995, but for other equiva- 
lent treatments, see, e.g., Bhattacharya & van den Heuvel 
1991; Nelemans et al. 1999). The eccentricity of the post- 
supernova binary is given by 
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where AMsn = ^He — ^bh' post-supernova semi- 
major axis by 

ao 
1 — e 

where ao is the initial orbital separation. The post- 
supernova system kick velocity can be obtained from equa- 
tion (2.10) in Brandt & Podsiadlowski (1995) as 
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where v^^^ is the pre-supernova orbital velocity of the sys- 
tem. Here we have neglected the small change in the mass 
of the secondary due to the capture of ejected material 
from the primary (typically ^0.2 M©), as well as any kick 
associated with the interaction of the supernova blast wave 
with the secondary (see Marietta et al. 2000). 
After the supernova, the binary parameters will continue 
to evolve. The system will first re-circularize, obtaining a 
new orbital separation 

arec = CipSN (1 " e^). 

Once the secondary starts to fill its Roche lobe, it will start 
to lose mass, of which a fraction /3 will be accreted by the 
primary, while the rest will be ejected from the system. We 
assume that the matter that is lost from the system carries 
away the same specific angular momentum as the primary 
(see, e.g., Podsiadlowski Rappaport, & Pfahl 2001). This 
is appropriate if the mass loss occurs from a region near 
the primary, as suggested by the relativist ic jets observed 
from Nova Sco (Hjellming & Rupen 1995). 
Even though this model is still relatively simple (for ex- 
ample, it does not take into account a kick due to an 
asymmetric explosion), it still contains a large number of 
essentially unspecified parameters (M^^, Mgj^, Mfaiiback, 

/ejection, /fallback, ^, P)- For givCU valuCS of /ejection and 

/fallback, wc havc Sampled all the other parameters in a 
fairly systematic and comprehensive fashion, although we 
generally do not change the present masses of the Nova Sco 
components, but keep them fixed at 5.4 and 1.45 M©, 
respectively^. In practice, we proceed in the following way. 
For each of the 7 supernova models (i.e., each combination 
of helium star mass and explosion energy), we systemat- 
ically vary the initial black- hole mass, Mgj^, the fallback 
mass, Mfaiiback, the wind- loss parameter, and the mass- 
accretion parameters, (3 (the latter two are varied from 
to 1). Having fixed these parameters, we can use the 
present orbital period and masses to reconstruct the pre- 
supernova masses and pre-supernova orbital period using 
the formalism outlined above. If this reconstruction shows 
that the radius of the pre-supernova secondary is smaller 
than its Roche lobe and that the system remains bound 
in the supernova explosion (if e < 1), we calculate the 
pollution and mixing in the secondary for our assumed 
values of /ejection and /fallback- In order to decide whether 
an acceptable model has been found, we define a quality 
parameter as 

2 



1=1 



[Xj/R] - [Xj/R] 

1 obs 
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where [X^/H] are the calculated logarithmic abundances 
(relative to solar) of N, O, Mg, Si, S, Ti, and Fe and 



^ We have also performed some calculations using masses of 6 Mq and 2Mq, respectively, obtaining similar results. 
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[X^/H]° ® are the abundances obtained by Israelian et al. 

(1999) and A[Xi/H]°^^ are the observational errors (see 
Table 1). We consider a particular model acceptable if 
Q < 1. 

2.2. Results 

In Table 2a we present the results for three combinations 
of /ejection and /fallback ([1,1], [0,1], [0,2]), for ah super- 
nova models that produced acceptable fits. The first two 
columns specify the initial mass of the helium star and 
the supernova explosion energy, while the next 10 columns 
give the mean logarithmic abundances (relative to solar) 
of He, C, O, Ne, Mg, Si, S, Ca, Ti, and Fe for acceptable 
models. The next 6 columns contain the mean values of 
selected model parameters (the mass of the He star just 
before the supernova, the initial mass of the compact star, 
the black-hole mass after fallback, the initial mass of the 
secondary, the system kick velocity and the maximum sys- 
tem kick velocity for each set of calculations, i.e., each of 
the 7 supernova models). The last figure A^tot gives the 
total number of acceptable models for each supernova set 
and provides a measure of how easy it is to obtain an 
acceptable model for the pollution in the secondary for 
each set. All the quoted uncertainties are the standard 
deviations calculated for all acceptable models in each set. 
They are not proper statistical error estimates, since they 
are based on an even sampling of the unknown parame- 
ters. Nevertheless they give an indication of the range of 
the corresponding model parameters. 
Consistency with nucleosynthesis calculations 
As the table shows, many acceptable fits can be obtained 
for many plausible combinations of the parameters. In 
particular, acceptable fits can be found for He star models 
with of 10 and 16 Mq for all explosion energies. To 
some degree, this just reflects the fact that explosive nu- 
cleosynthesis produces similar overall abundance patterns 
in all of these models. Indeed, the overall consistency of 
the modelled abundances with the observed pollution in 
the secondary of Nova Sco provides confirmation of the 
general predictions of stellar and supernova nucleosynthe- 
sis. 

The only exception to this picture is the abundance of Ti, 
which is too low by at least a factor of 2 in all models. 
This may be caused by errors in some of the nuclear cross 
sections used in the explosive nucleosynthesis calculation 
or could provide evidence for a more a complicated nucle- 
osynthesis environment in the explosion (see § 3.5). 
In Table 3 we summarize some of the key binary param- 
eters of the best-fit models for Models A, B and C in 
Table 2a. Note particularly that the initial separation 
before the supernova is typically less than 6 i?©, much 
smaller than the present separation (15.2 i?©). This im- 
plies that the secondary was almost filling its Roche lobe 
before the supernova and results in a relatively large geo- 
metrical cross section for the capture of material from the 
supernova ejecta (^ 0.26 Mq). The semi-major axis of the 
post-supernova system is only slightly increased (at least 
in models where the black hole receives no kick at birth). 
Hence there had to be substantial mass transfer (^ 1 Mq) 
after the supernova to widen the system to the present 
separation. 

The model also predicts some other abundance anomalies. 



The abundance of C is predicted to be enhanced by up to 
a factor of 1.7, the abundance of Ca by a factor of 3 to 7, 
and the abundance of Ne by a factor of 2 to 5. The en- 
hancement of Ca is weakly correlated and the enhancement 
of Ne is weakly anti-correlated with the explosion energy 
(these correlations are strongest for the 10 Mq models). 
Therefore, these patterns provide, at least in principle, a 
means by which one could distinguish between a normal 
supernova and a hypernova event. 
Requirement of fallback 

All models require some fallback, where hypernova models 
require the smallest amount. This is the result of two fac- 
tors: (1) the initial black-hole mass must be close to the 
mass cut since matter near the mass cut (which is much 
smaller than the present black-hole mass) must have been 
ejected in order to produce the observed abundances. The 
mass cut is a function of the mass of the helium star and 
the energy of the explosion and is largest for the 16 
hypernova model (see Fig. 1). (2) The black- hole mass 
cannot increase by much more than ^ 1 Mq by mass ac- 
cretion from the secondary after the supernova. This is 
a consequence of the constraints on the pre-supernova or- 
bital parameters imposed by the present orbital separa- 
tion and component masses. A larger amount of post- 
supernova mass transfer implies a tighter pre-supernova 
binary (because mass transfer widens the system). The 
largest amount of mass transfer is therefore determined 
by the condition that the pre-supernova secondary just 
fills its Roche lobe. We note that this constraint would be 
substantially weaker if we had allowed for an asymmetric 
supernova explosion. 
Low system kick velocities 

The typical system kick velocities are relatively low, vary- 
ing from 10 to 60kms~^. There are a few extreme cases 
for the 16 Mq models where the system velocity is as high 
as ~ 90kms~^. However, these cases are very rare and 
require very special model parameters. These are cases 
where the helium star loses relatively little of its envelope 
in a stellar wind before the supernova (Mne ^13 Mq) and 
where the system becomes almost unbound in the super- 
nova event (the immediate post-supernova eccentricity is 
close to 1). Whether such a small amount of wind mass loss 
is physically reasonable is somewhat doubtful. In any case, 
since the parameter range that leads to these high kick ve- 
locities is extremely limited, such solutions, while not for- 
mally ruled out, are statistically not favoured. The min- 
imum observed space velocity of Nova Sco is 106kms~^ 
and possibly much larger, since this corresponds to the ra- 
dial velocity only. It therefore appears very unlikely that a 
symmetric supernova explosion alone could explain the ob- 
served velocity. This implies that an additional kick, e.g., 
due to an asymmetry in the explosion, is required. While 
our model did not take this possibility into account, our re- 
sults are consistent with it. The initial masses of the com- 
pact remnant, Mgj^, for many of the helium-star models 
with a mass of 10 Mq and for normal supernova explosions 
are consistent with the maximum mass allowed for neutron 
stars. Even when M^^ is larger than the maximum mass 
of a neutron star, a neutron star may have formed first, 
since, in our simple model, Mgj^, strictly speaking, does 
not only include the initial mass of the compact remnant, 
but also any fallback material that did not reach the or- 
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bit of the secondary. Therefore, a two-stage process where 
the cohapse first leads to the forraation of a neutron star, 
accompanied by a supernova kick, which is subsequently 
converted into a black hole by fallback may well explain 
the observed space velocity (see Brandt et al. 1995). 
We note that these conclusions are not inconsistent with 
the findings of Nelemans et al. (1999). Their results also 
imply that one requires rather special parameters to ex- 
plain the observed space velocity of Nova Sco with a sym- 
metric supernova explosion alone. In addition, some of the 
sets of parameters they present for illustration can be ruled 
out by the present investigation which takes into account 
more observational constraints and the reduced masses of 
the Nova Sco components. 

3. DISCUSSION 

In the previous section we showed that a relatively sim- 
ple model can explain most of the observed a-element en- 
hancements in the secondary of Nova Sco. We now turn 
to an examination of the question whether some of the as- 
sumptions that went into the model are actually physically 
reasonable. 

3.1. Capture efficiency 

Obviously, a significant fraction of the supernova ejecta 
must have been captured by the secondary. This is not 
necessarily expected. Marietta et al. (2000) recently ex- 
amined the effects of a supernova blast wave on a binary 
companion in the context of a Type la supernova. They 
showed that a significant part of the outer part of the 
secondary was lost by ram-pressure stripping. Instead of 
accreting matter from the supernova ejecta, the secondary 
actually lost mass. On the other hand, the layers in which 
a-process elements are synthesized are buried deep inside 
the core and move with a much lower velocity than the 
outer layers of the helium star, which are mainly respon- 
sible for the ram-pressure stripping. This makes it easier 
for some of that material to be captured by the secondary 
(see the discussion in Marietta et al. 2000). 
The problem is less severe if the material captured comes 
mainly from fallback since this material will generally have 
a much lower velocity than the escaping material. Indeed, 
if the velocity relative to the secondary is smaller than the 
escape velocity at the surface of the secondary, gravita- 
tional focusing could increase the capture efficiency above 
a value of 1. 

3.2. The requirement of fallback 

Substantial fallback is commonly found in supernova ex- 
plosions with hydrogen-rich envelopes (see Woosley & 
Weaver 1995). This fallback occurs either because matter 
(typically just above the mass cut) did not attain escape 
speed in the supernova or, more significantly, because it 
was pushed back by hydro dynamical effects, in particu- 
lar by a reverse shock that forms at the interface of the 
He/H boundary. However, in the present model there is 
no hydrogen envelope to produce the interface for a strong 
reverse shock. A strong reverse shock is also formed at 
the CO/He interface as seen in Type lb supernova models 
(Hachisu et al. 1991, 1994). In the absence of a strong 
reverse shock from the composition interface, the natural 
location for fallback, if it occurs, is close to the original core 



and not close to the location of the secondary (as required 
in our simple model) (unless the system is still surrounded 
by a common envelope, which cannot be completely ruled 
out, but is a priori very unlikely). 

3.3. Fallback with mixing? 

A possible solution to this problem is a modification of 
our model where fallback occurs from a region near the 
core, but is accompanied by substantial mixing. Some of 
this mixed material, enriched with a-element-rich material 
near the mass cut, is ejected and produces the observed 
pollution in the secondary, while the rest falls back onto 
the compact remnant (see Fig. 2, Case B). 
There are several scenarios in which such mixing could 
arise. (1) In some of the collapsar model studied by Mac- 
Fadyen, Woosley & Heger (2001), the initial collapse leads 
to the formation of a neutron star and a weak supernova 
shock which fails to eject all of the helium core. As the 
shock stalls, matter falls back onto the core converting 
the neutron star into a black hole. The energy that is 
associated with this collapse may then be able to eject 
the outer envelope (perhaps in the form of jet-powered 
shocks; for a similar suggestion for Nova Sco, see Brown 
et al. 2000). Vigorous mixing is expected to occur in this 
scenario, since, as material tries to fall back, it is be- 
ing pushed up by a low-density neutrino-heated bubble, 
which is Ray leigh- Taylor unstable (see, e.g., Kifonidis et 
al. 2000). 

(2) If a black hole forms promptly, mixing in the ejecta 
may be induced by Ray leigh- Taylor instabilities at the 
Si/0 interface (Kifonidis et al. 2000; Kifonidis 2001) as 
weh as the CO/He interface (Hachisu et al. 1991, 1994). 
After mixing, the inner part of the mixed layers may fall 
back, increasing the mass of the initial black hole. Umeda 
& Nomoto (2001) have shown that such mixing followed 
by fallback can explain the large Zn abundance observed 
in very metal-poor stars and suggested that this may be a 
generic feature of core-collapse supernovae, either involv- 
ing the formation of a black hole or a neutron star. 
3) In an aspherical explosion, heavy elements synthesized 
in a deep layer may be mixed into outer layers in the form 
of a jet (Maeda et al. 2001), while fallback occurs from the 
equatorial region. This is equivalent to the mixing & fall- 
back process discussed in Umeda & Nomoto (2001) and 
could occur for both the formation of a black hole or a 
neutron star. 

Such a scenario has several advantages: (1) a neutron star 
may form first, and the associated neutrino emission can 
be responsible for a standard neutrino-induced supernova 
kick to explain the observed space velocity. (2) Material 
near the mass cut is mixed into the envelope and can es- 
cape with it, explaining the abundance anomalies in the 
secondary. (3) The supernova explosion could be weaker 
and ram-pressure stripping would then be less of a prob- 
lem; this could allow for more efficient capture of super- 
nova material by the secondary (even in the case of a hy- 
pernova explosion, the ejecta in the equatorial plane could 
be moving relatively slowly). 

3.4. Case B: Fallback with mixing simulations 

In order to simulate a fallback model with mixing, we mod- 
ified the model described in § 2.1 by assuming that the 
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layer between Mgj^ and Mgjj+AMmix is completely mixed 
during the collapse phase, where we defined AMmix = 
^ ^fallback (see Fig. 3b). The definition of includes 
both the initial mass of the compact remnant and any 
fallback material that was not mixed with the material 
that is ejected (this could, for example, be material from 
the inner parts of the disk that forms around the initial 
compact core, as in the collapsar models of MacFadyen & 
Woosley (1999) and MacFadyen et al. (2001) or material 
falling back from the equatorial region of the star in an 
aspherical explosion (Maeda et al. 2001). 
In Table 2b we present the results of these simulations for 
m = 1, 2, 3, and 4. Most of the results are similar to 
the earlier calculations: the amount of fallback is similar, 
while the system velocities tend to be a bit larger. One 
significant difference is that, with the inclusion of mixing, 
hypernova models are preferred over standard supernova 
models, which tend not to produce enough S. 
The size of the mixing region outside the initial compact 
remnant varies quite significantly between different mod- 
els; in some hypernova models, the mixing region contains 
less than 1 Mq , while in normal supernova models at least 
several Mq are required. 

3.5. Case C: Non- spherical models 

So far, we have assumed that the ejection of matter occurs 
in a more-or-less spherically symmetric way. However, this 
is not generally expected for hypernova or collapsar mod- 
els, in particular those associated with gamma-ray bursts 
(MacFadyen & Woosley 1999). In these models it is gen- 
erally believed that the core is initially rapidly rotating 
and that the accretion of matter onto the compact ob- 
ject occurs through an accretion disk. Maeda et al. (2001) 
have recently constructed aspherical hypernova models for 
the hypernova SN 1998bw, simulating both the hydrody- 
namics and the nucleosynthesis in two dimensions for a 
helium- star model of 16 Mq at the beginning of helium 
burning. They showed that the chemical composition of 
the ejecta is strongly dependent on direction. In particu- 
lar, Fe is mainly ejected in the polar direction, while O and 
Mg are preferentially ejected near the equatorial plane. 
To simulate the pollution by an aspherical hypernova, we 
assume that the secondary is located in the equatorial 
plane of the helium star (and the black- hole accretion disk) 
and that the secondary captures material that is within an 
angle 6 of the equatorial plane, where 6 is the angular ra- 
dius subtended by the secondary as seen from the helium 
star. The results of these simulations (for their best model 
C) are presented in Tables 3 and 4 (models indicated by 
a *). Somewhat surprisingly, none of our simulations pro- 
duced acceptable fits (the results shown were obtained by 
increasing the acceptance parameter Q from 1 to 2). This 
is a direct consequence of the large overabundance of O 
and Mg near the equatorial plane. All models produce ei- 
ther an unacceptable overabundance of O and Mg or an 
unacceptable underabundance of S and Si, depending on 
where the cut-off below which matter can be mixed into 
the ejecta occurs. 

However, when modeling the shapes of spectral lines in SN 
1998bw, Maeda et al. (2001) found that the fits could be 
improved if there was some lateral mixing in the ejecta, 
i.e., between material ejected in the equatorial plane and 



material ejected more along the jet axis (e.g., due to a 
shear instability). To test this possibility we also con- 
sidered a model where we assumed complete lateral mix- 
ing in the ejecta for the same hypernova model as used 
above, i.e., where all the material within given velocity 
bins was assumed to be mixed completely. The results 
of these simulations are also shown in Tables 3 and 4 
and indicated with a prime ('). In this case, excellent 
fits are obtained, as one would have expected, since this 
model should approximate the spherical hypernova model 
used earlier. Interestingly, the Ti abundance is signifi- 
cantly enhanced ([Ti/H]~ 0.5) and is marginally consis- 
tent with the observed value ([Ti/H]- 0.9 =b 0.4). This 
increased Ti abundance comes at the price of an increased 
Fe abundance. While both the Ti and Fe abundances are 
consistent with the observed values, the ratio [Ti/Fe] is 
still significantly below the observed one (^ 0.2 instead 
of 0.8). We emphasize that the assumption of complete 
lateral mixing is extreme; at present we can not identify 
any physical process that would lead to such a result. We 
note, however, that Hofiich, Khokhlov & Wang (2001) have 
shown that the jet in an aspherical supernova model is de- 
celerated at the H/He interface and that material spreads 
laterally, although only to a limited extent. Whether this 
could provide a viable model, requires computations with 
much higher numerical resolution. Irrespectively, our re- 
sults suggest that the enhanced abundance of Ti (which 
results from a mixture of the nucleosynthesis products of 
complete and incomplete Si burning; see Umeda & Nomoto 
[2001] for a discussion) could potentially provide a signa- 
ture for an asymmetric hypernova. 

4. CONCLUSIONS 

The main conclusion of this investigation is that, us- 
ing standard supernova models and a relatively simple 
model for the pollution of the secondary, we can explain 
the observed a-element enhancements in the secondary of 
Nova Sco, confirming standard nucleosynthesis predictions 
(apart from the abundance of Ti, which is always too low). 
Nova Sco presents a clear case for a two-step black-hole for- 
mation process, where a substantial fraction of the black- 
hole mass is the result of fallback. In order for the sec- 
ondary to be polluted with material near the mass cut, 
this fallback material must either have reached the loca- 
tion of the secondary before falling back or, more likely, 
be mixed during the explosion with material that will es- 
cape. A two-step black-hole formation process may pro- 
vide a simple explanation for the high space velocity of 
Nova Sco, since the system may have received the same 
type of neutrino-induced kick as normal neutron stars are 
believed to. Alternatively, asymmetric mass ejection in 
an aspherical super nova/hyper nova model may also con- 
tribute to the observed space velocity. 
Our analysis shows that helium star models of 10 to 16 Mq 
are most probable and that both normal supernova as well 
as more energetic hypernova models can explain the abun- 
dances observed in the secondary. The majority of accept- 
able models are, however, hypernova models, in particular 
for the more realistic models that include mixing. The 
one aspherical hypernova model we used only provided an 
acceptable fit to the observed abundance anomalies when 
we assumed that there was extreme lateral mixing between 
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the eject a in the equatorial plane and in the polar direc- 
tion. Of course, so far we only tried one helium-star model 
with a mass of 16 Mq for one particular explosion energy. 
It is quite possible that, for a lower-mass model, the oxy- 
gen produced in the equatorial plane would be lower, while 
still preserving most of the hypernova features without re- 
quiring extreme lateral mixing. The modeling of the pol- 
lution for the aspherical explosion is clearly even more un- 
certain than in the spherical case, which suggests that it is 
important to study mixing processes in aspherical models 
with much higher numerical resolution. 
On the other hand, it is also worth noting that it is not 
necessarily expected that the helium-star progenitor of the 
black hole was rapidly rotating, since the progenitor is 
likely to have passed through an extended red-supergiant 
phase before the system experienced a common-envelope 
and spiral-in phase. It is quite likely that, in the red- 
supergiant phase, the helium core would have been sig- 
nificantly spun down by both hydrodynamical (Heger, 
Langer, & Woosley 2000) and magnet ohydro dynamical 
processes (Spruit & Phinney 1998), which efficiently cou- 
ple the core to the slowly rotating, convective envelope 
(as may also be required to explain the low initial spin 
periods of the majority of radio pulsars in supernova rem- 
nants). Strohmayer (2001) has provided some evidence 
that the black hole may be rotating rapidly at the present 
time (based on observations of quasi-periodic oscillations 
at 450 Hz). However, as our modeling has shown, the black 
hole is likely to have accreted ^ 1 Mq after the supernova 
from the companion star (also see Beer & Podsiadlowski 



2001). Since this accretion would also have spun up the 
black hole, the observations by Strohmayer (2001) do not 
provide a strong constraint on the rotation of the black 
hole at birth. 

It is interesting to speculate on the relation between Nova 
Sco and other low-mass black-hole binaries. Most low- 
mass black- hole binaries have rather low space velocities, 
entirely consistent with the expected velocity dispersion 
caused by scattering by molecular clouds and spiral arms 
(Brandt et al. 1995). This may suggest that the black 
holes in these systems formed promptly, i.e., without a 
kick or significant mass ejection (see, however, the discus- 
sion in Nelemans et al. 1999 for a rather different point 
of view). This would then imply that the secondaries in 
these systems should not show the same enhancement in 
ce-process elements as the secondary in Nova Sco, a pre- 
diction that should be checked with future high-resolution 
spectral observations of these systems. Orosz et al. (2001) 
recently reported the discovery of similar over-abundances 
of a-process elements in the B-star secondary of another X- 
ray transient J1819.3-2525 (V4641 Sgr). Interestingly, this 
systems also appears to have an anomalously high space 
velocity (its measured 7- velocity is 107kms~^). Unfortu- 
nately, this velocity estimate is quite uncertain at present 
because of the proximity of the system to the Galactic cen- 
ter, which makes the distance estimate very sensitive to 
the precise distance. Nevertheless, if confirmed, it would 
suggest that the system may have experienced a similar 
evolutionary history as Nova Sco. 
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1 . Wide binary with massive 
primary star and intermediate- 
mass companion. 



2. Spiral-in of secondary in red supergiant 
phase leads to ejection of the hydrogen-rich 
envelope and produces short-period helium 
star binary. 



Case 




3ab. Core collapse 
creates neutron star 
and leads to ejection 
of primary's envelope, 



4b. Fallback occurs in region 
interior to secondary's orbit, with 
significant mixing of the material 
at the fallback interface. 
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Case C 




::::0 



3c. Asymmetric core 
collapse creates neutron 
star/black hole and initiates 
aspherical expansion of the 
primary's 



4c. Mixing 
occurs on the 
fallback 
interface. The 
material in 
the 

equatorial 
region 
expands 
more slowly 
than along the 
poles. 




5c. Fallback 
increases mass 
of compact 
object. Slow- 
moving 
equatorial 
debris 
deposits 
synthesized 
material on 
secondary. 



Fig. 2. — Schematic diagram for the various pollution scenarios considered. 
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Fig. 3. — Schematic diagram defining the various mass parameters in the simple fallback model (a) and the fallback with mixing model (b). 
M£pj is the initial black-hole mass, the black-hole mass after fallback, Mne the mass of the helium star at the time of the supernova, 

-^He initial mass of the helium core (without wind mass loss), AMmix the mass in the mixing region (b). 
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TABLE 1 

Observed abundances in the secondary of Nova Sco^ 





N 





Mg 


Si 


s 


Ti 


Fe 


A [X,/H] 


0.45 
0.50 


1.00 
0.30 


0.90 
0.40 


0.90 
0.30 


0.75 
0.20 


0.90 
0.40 


0.10 
0.20 



Note. — [X^/H]: logarithmic abundances relative to solar; 
A[X^/H]: observational uncertainties in [X^/H]. 
References. — (1) Israelian ct al. 1999 



TABLE 2a 

Acceptable fallback models for the secondary of Nova Sco 



(Me) (foe) 


He 


C 





Ne 


Mg 


Si 


S 


Ca 


Ti 


Fe 


(^0) 


(^0) 


(M0) 


(^0) 


^kick Cck 

(km/s) (km/s) 


A^tot 


Model A: Capture models with 


fallback 


(/ejection 


— 1, /fallback — 


1, M2 = 


1.45 Mq 


,Mbh 


= 5.4 Mq) 














1 n 


1 


0.03 


0.21 


0.76 


0.69 


0.83 


0.84 


0.73 


0.53 


0.22 


0.20 


6.96 


2.04 


4.66 


2.36 


39 


81 


79242 






ibO.Ol 


±0.03 


±0.06 


±0.06 


±0.06 


±0.09 


±0.09 


±0.10 


±0.06 


±0.14 


±1.74 


±0.18 


±0.50 


±0.45 


it io 






16 


1 


0.02 


0.11 


0.92 


0.65 


0.89 


0.85 


0.74 


0.59 


0.19 


0.19 


8.21 


2.74 


4.93 


2.04 


46 


90 


64780 






lb 0.00 


±0.03 


±0.10 


±0.09 


±0.10 


±0.10 


±0.11 


±0.14 


±0.08 


±0.14 


±2.03 


±0.23 


±0.34 


±0.35 


it ly 






10 


8 


0.03 


0.20 


0.72 


0.53 


0.73 


0.88 


0.81 


0.69 


0.17 


0.15 


6.95 


2.63 


4.65 


2.38 


40 


81 


61514 






±0.00 


±0.03 


±0.06 


±0.05 


±0.06 


±0.08 


±0.10 


±0.14 


±0.07 


±0.10 


±1.75 


±0.16 


±0.51 


±0.45 


it io 






10 


30 


0.03 


0.21 


0.65 


0.36 


0.62 


0.91 


0.85 


0.74 


0.15 


0.10 


7.10 


3.61 


4.68 


2.45 


44 


81 


16952 






±0.00 


±0.02 


±0.03 


±0.02 


±0.02 


±0.06 


±0.09 


±0.13 


±0.05 


±0.05 


±1.74 


±0.14 


±0.49 


±0.40 


J- on 
it zU 






16 


30 


0.02 


0.09 


0.81 


0.51 


0.74 


0.89 


0.84 


0.80 


0.20 


0.16 


8.42 


4.22 


4.99 


2.00 


47 


90 


69290 






±0.00 


±0.02 


±0.10 


±0.08 


±0.10 


±0.10 


±0.12 


±0.16 


±0.08 


±0.09 


±1.99 


±0.27 


±0.28 


±0.31 


±19 






Model 


B: Capture models with 


fallback 


( t/ccLpture 


— 0, /fallback — 


1, M2 = 


1.45 M® 


, Mbh ■ 


= 5.4 Mq) 














10 


1 


0.03 


0.24 


0.68 


0.62 


0.73 


0.75 


0.64 


0.47 


0.20 


0.22 


4.48 


1.92 


3.74 


3.27 


22 


45 


6996 






±0.00 


±0.01 


±0.02 


±0.03 


±0.03 


±0.02 


±0.02 


±0.05 


±0.05 


±0.10 


±0.69 


±0.12 


±0.62 


±0.40 


±6 






16 


1 


0.02 


0.12 


0.87 


0.63 


0.85 


0.77 


0.65 


0.54 


0.19 


0.24 


6.95 


2.58 


4.66 


2.38 


42 


90 


28655 






±0.00 


±0.02 


±0.08 


±0.08 


±0.09 


±0.05 


±0.06 


±0.07 


±0.08 


±0.14 


±1.49 


±0.21 


±0.37 


±0.29 


±17 






10 


8 


0.03 


0.24 


0.64 


0.53 


0.66 


0.76 


0.71 


0.66 


0.18 


0.19 


4.48 


2.48 


3.79 


3.25 


20 


36 


4200 






±0.00 


±0.01 


±0.02 


±0.05 


±0.03 


±0.03 


±0.03 


±0.04 


±0.03 


±0.07 


±0.65 


±0.08 


±0.63 


±0.38 


±5 






10 


30 


0.03 


0.23 


0.59 


0.38 


0.61 


0.81 


0.77 


0.72 


0.29 


0.24 


4.40 


3.01 


3.81 


3.34 


18 


26 


1317 






±0.00 


±0.01 


±0.03 


±0.01 


±0.02 


±0.05 


±0.08 


±0.10 


±0.14 


±0.10 


±0.72 


±0.35 


±0.76 


±0.48 


±7 






16 


30 


0.02 


0.10 


0.77 


0.55 


0.76 


0.76 


0.71 


0.74 


0.24 


0.26 


7.00 


3.85 


4.82 


2.29 


38 


84 


32257 






±0.00 


±0.02 


±0.09 


±0.06 


±0.08 


±0.06 


±0.08 


±0.10 


±0.10 


±0.13 


±1.29 


±0.31 


±0.27 


±0.23 


±17 







TABLE 2a— Continued 



ML 




He 


C 





Ne 


Mg 


Si 




Ua 


il 


Fe 






/\/f 1 
^^BH 


M2 


^kick ^ 


jinax 
kick 


AT 


(M@) (foe) 
























(^©) 


(^©) 


(^©) 


(km/s) (k: 


m/s) 




Model C: Capture models with fallback (/capture 


— 0, /fallback — 


I, M2 


1 /I 1\ /f 

1.45 Mq 




= 5.4 M®) 














10 


1 


0.02 


0.16 


0.69 


0.61 


0.76 


0.85 


0.75 


0.54 


0.21 


0.17 


6.58 


2.07 


4.57 


2.48 


37 


81 


54039 






ibO.Ol 


±0.07 


±0.06 


±0.06 


±0.06 


±0.08 




-1- n 1 n 




±0.12 


±1.74 




ziz yj.o^ 


ZIZ yj.^o 


-1- 1 Q 
ZIZ ±y 






16 


1 


0.02 


0.10 


0.84 


0.57 


0.81 


0.85 


0.76 


0.61 


0.19 


0.18 


7.73 


2.77 


4.85 


2.14 


44 


90 


49459 






lb 0.00 


±0.02 


±0.11 


±0.10 


±0.11 


±0.09 


±0.11 


±0.14 


±0.08 


±0.14 


±1.87 


±0.21 


±0.34 


±0.32 


±19 






10 


8 


0.03 


0.19 


0.67 


0.45 


0.66 


0.90 


0.84 


0.73 


0.17 


0.15 


5.78 


2.63 


4.44 


2.69 


27 


81 


22673 






ibO.Ol 


±0.06 


±0.05 


±0.08 


±0.05 


±0.07 


±0.08 


±0.12 


±0.07 


±0.08 


±1.48 


±0.13 


±0.62 


±0.45 


±15 






10 


30 


0.03 


0.22 


0.61 


0.38 


0.62 


0.83 


0.78 


0.72 


0.22 


0.18 


4.51 


3.37 


3.92 


3.18 


17 


26 


1617 






lb 0.00 


±0.01 


±0.03 


±0.02 


±0.02 


±0.07 


±0.08 


±0.10 


±0.08 


±0.10 


±0.67 


±0.27 


±0.69 


±0.35 


±5 






16 


30 


0.02 


0.10 


0.75 


0.53 


0.73 


0.85 


0.82 


0.83 


0.21 


0.19 


6.99 


4.16 


4.82 


2.23 


37 


81 


30943 






lb 0.00 


±0.02 


±0.10 


±0.07 


±0.10 


±0.09 


±0.12 


±0.13 


±0.08 


±0.09 


±1.24 


±0.22 


±0.29 


±0.25 


±16 







Note. — Mpj^: initial helium-star mass; E^n'- explosion energy (Ifoe = 10^^ erg); He, C,. . . Fe: logarithmic abundances (relative to solar) of He, C,. . . Fe; Mhc^ 
pre-supernova helium-star mass; M^^: initial mass of the compact remnant; M^^: black-hole mass after fallback; : initial mass of the secondary; t^kick- system 
kick velocity; v^^^: maximum system kick velocity for each supernova set; A^tot- total number of acceptable models for each supernova set. 



TABLE 2b 

Acceptable mixing models for the secondary of Nova Sco 



AfO, EsN He C O Ne Mg Si S Ca Ti Fe M^e AM^i, v^i,^ v^^^ Nu 

(M@)(foe) (M@) (M@) (M@) (M@) (M@) (km/s) (km/s) 



Model 
10 

16 

Model 
16 



D; Capture models with mixing 

30 0.03 0.23 0.63 

±0.00 ±0.00 ±0.01 

30 0.02 0.12 0.93 

±0.00 ±0.01 ±0.05 

E: Capture models with mixing 



1 



10 
16 



30 



30 



0.15 

±0.01 

0.23 

±0.00 

0.11 

±0.01 



0.97 

±0.03 

0.63 

±0.01 

0.86 

±0.07 



0.03 
±0.00 
0.03 
±0.00 
0.02 
±0.00 

Model F: Capture models with mixing 

0.88 
±0.02 

0.69 
±0.00 
0.63 
±0.01 
0.84 
±0.07 



16 


1 


0.02 


0.12 






±0.00 


±0.01 


10 


8 


0.03 


0.24 






±0.00 


±0.00 


10 


30 


0.03 


0.23 






±0.00 


±0.00 


16 


30 


0.02 


0.11 






±0.00 


±0.02 



; (/capture 

0.39 
±0.00 

0.62 
±0.04 

(/capture 

0.74 
±0.03 

0.39 
±0.00 

0.58 
±0.05 

(/capture 

0.65 
±0.02 

0.57 

±0.00 

0.39 

±0.01 

0.56 

±0.06 



Model G: Capture models with mixing (/captur< 

0.67 
±0.00 
0.39 
±0.01 
0.56 
±0.06 



16 


1 


0.02 


0.13 


0.90 






±0.00 


±0.00 


±0.00 


10 


30 


0.03 


0.23 


0.63 






±0.00 


±0.00 


±0.01 


16 


30 


0.02 


0.11 


0.85 






±0.00 


±0.02 


±0.07 



= 1, AM^i, = 
0.64 0.87 
±0.01 ±0.04 
0.87 0.83 
±0.04 ±0.09 



1, AM^ 



0.96 

±0.03 

0.63 

±0.01 

0.82 

±0.06 



0.70 
±0.02 

0.83 
±0.05 
0.75 
±0.07 



= 1, AM,,ix = 
0.87 0.69 



±0.02 

0.72 

±0.00 

0.63 

±0.01 

0.80 

±0.07 



±0.02 

0.70 

±0.00 

0.84 

±0.05 

0.77 

±0.07 



1, AM„ 



0.89 

±0.00 

0.63 

±0.01 

0.80 

±0.07 



0.70 

±0.01 

0.85 

±0.04 

0.79 

±0.07 



MMiback, M2 = 1 

0.83 0.77 

±0.06 ±0.08 

0.71 0.61 

±0.13 ±0.17 

2 Mfallback, M2 = 

0.50 0.42 

±0.01 ±0.02 

0.79 0.73 

±0.06 ±0.07 

0.63 0.61 

±0.07 ±0.08 

3MMiback, M2 = 

0.54 0.45 

±0.02 ±0.03 

0.61 0.52 

±0.01 ±0.02 

0.80 0.73 

±0.06 ±0.08 

0.68 0.66 

±0.07 ±0.09 

4 Mfallback, M2 = 

0.53 0.38 

±0.01 ±0.01 

0.81 0.75 

±0.06 ±0.08 

0.70 0.66 

±0.08 ±0.11 



.45 M®, Mbh = 

0.20 0.17 

±0.06 ±0.07 

0.13 0.10 

±0.07 ±0.07 

1.45 Mq, Mbh ^ 

0.15 0.19 

±0.03 ±0.07 

0.26 0.27 

±0.10 ±0.11 

0.31 0.28 

±0.18 ±0.14 

1.45 M®, Mbh ^ 

0.16 0.20 

±0.04 ±0.09 

0.12 0.10 

±0.01 ±0.01 

0.20 0.19 

±0.05 ±0.08 

0.27 0.26 

±0.17 ±0.15 

1.45 M®, Mbh 

0.08 0.06 

±0.00 ±0.00 

0.20 0.18 

±0.05 ±0.08 

0.19 0.19 

±0.10 ±0.13 



5.4 M®) 
3.94 
±0.15 
6.05 
±0.91 

= 5.4 M^ 
7.01 
±1.20 
3.98 
±0.28 
8.24 
±1.97 

= 5.4 Mc, 

10.71 
±0.84 

3.59 
±0.02 

4.13 
±0.37 
8.97 
±1.97 

= 5.4 Mf 

11.60 

±0.00 

4.04 

±0.24 

9.11 

±2.04 



2.57 


3.27 


3.58 


0.71 


22 


25 


7396 


-1-0/19 


-1- n 1 ^ 


-1- n 1 ^ 

It U.iO 


It U.^iZ 


It z 






3.40 


4.34 


2.59 


0.94 


37 


66 


27814 


It U.O / 


It U.Z4 


It U.zz 


It U.O / 


It io 






)) 

2.54 


4.43 


2.50 


3.78 


53 


76 


1828 


±0.14 


±0.35 


±0.28 


±0.84 


± 11 






2.91 


3.28 


3.57 


0.74 


23 


35 


4035 


±0.29 


±0.21 


±0.21 


±0.40 


±3 






3.44 


4.80 


2.17 


2.73 


54 


90 


80486 


-1- n ^/i 

□I U.O'i 


It u.oo 


It u.oz 


It i.OO 


It io 






)) 

2.55 


5.01 


1.90 


7.38 


80 


90 


2353 


±0.16 


±0.14 


±0.09 


±0.69 


±5 






2.55 


2.88 


3.97 


0.98 


26 


26 


36 


±0.02 


±0.02 


±0.02 


±0.06 


±0 






3.19 


3.40 


3.45 


0.61 


23 


45 


1694 


±0.14 


±0.20 


±0.20 


±0.46 


±5 






3.68 


4.89 


2.06 


3.63 


60 


90 


62403 


±0.50 


±0.33 


±0.31 


±1.98 


±17 






2.84 


5.01 


1.84 


8.67 


90 


90 


11 


±0.02 


±0.00 


±0.00 


±0.06 


±0 






3.23 


3.35 


3.50 


0.48 


23 


35 


1052 


±0.14 


±0.16 


±0.16 


±0.32 


±3 






3.95 


4.89 


2.05 


3.74 


61 


90 


39602 


±0.36 


±0.33 


±0.32 


±2.13 


±17 







Note. — AM^ix: mass of the mixing region above Mgp^; all other columns are as in Table 2a. 
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TABLE 3 

Binary parameters for "best-fit" models 



Model 


ML 

{Mr,) 


(foe) 


Mue 

{Mr,) 


(Me) 


Mhn 

{Mr,) 


M^ 

{Mr,) 


{R&) 


'^post 

{R@) 


e 


/geo 


-^capture 


^kick 

(km/s) 


Fallback Models 
























A 


10.00 


1 


5.57 


2.22 


5.07 


2.76 


5.68 


6.06 


0.06 


0.025 


0.27 


11 


B 


16.00 


30 


5.56 


3.94 


4.40 


2.78 


5.84 


6.96 


0.16 


0.024 


0.26 


28 


C 


16.00 


30 


5.54 


4.06 


4.38 


2.81 


5.79 


6.91 


0.16 


0.025 


0.28 


28 


A* 


15.88 


10 


4.96 


2.40 


3. 75 


3.10 


5.84 


7.10 


0.18 


0.029 


0.51 


35 


B* 


15.88 


10 


5.50 


2.40 


4.35 


2.85 


5.73 


6.82 


0.16 


0.026 


0.37 


29 


A' 


16.00 


10 


6.20 


4.30 


5.11 


2.61 


5.66 


6.59 


0.14 


0.023 


0.29 


23 


B' 


16.00 


10 


6.61 


2.48 


4.26 


2.59 


5.68 


8.65 


0.34 


0.023 


0.25 


54 


Mixing 


; Models 
























D 


16.00 


30 


3.81 


3.80 


3.81 


3.04 


6.93 


6.93 


0.00 


0.020 


0.24 





E 


16.00 


30 


6.63 


3.81 


4.28 


2.57 


5.73 


8.70 


0.34 


0.022 


0.26 


53 


E 


16.00 


30 


6.89 


3.89 


4.61 


2.50 


5.66 


8.33 


0.32 


0.022 


0.25 


48 


G 


16.00 


30 


8.00 


3.94 


4.57 


2.28 


5.71 


11.43 


0.50 


0.018 


0.21 


65 




15.88 


10 


6.58 


2.40 


4.26 


2.60 


5.69 


8.61 


0.34 


0.023 


0.30 


53 




15.88 


10 


10.45 


2.40 


5.01 


1.96 


5.56 


25.08 


0.78 


0.015 


0.20 


80 


E' 


16.00 


10 


6.61 


2.46 


4.26 


2.59 


5.68 


8.65 


0.34 


0.023 


0.26 


54 




16.00 


10 


9.24 


2.48 


174 


2.12 


5.59 


16.33 


0.66 


0.017 


0.19 


76 



Note. — The present masses of the black hole and the secondary in the Nova Sco system are assumed to be 5.4 
and 1.45 M^g), respectively, implying a present orbital separation of 15.2 i^0. ttp^e/ttpost- pre-/post-supernova semi- 
major axis; e: post-supernova eccentricity; /geo- pre-supernova fractional, geometrical cross section of the secondary; 
Mcapture^ total mass captured by the secondary; all other parameters are the same as in Tables 2a and 2b. 



TABLE 4 

Best fallback/mixing models for aspherical hypernova models 



Mg, £'SN He C O Ne Mg Si S Ca Ti Fe Mne ^2 ^^mix %ick <eT ^t, 

(Me) (foe) (Me) (Me) (Me) (M@) (M©) (km/s) (km/s) 



A^o Lateral Mixing 

Model A"" : Capture models with fallback (/ejection 
16 10 0.02 0.16 1.12 0.79 0. 

±0.00 ±0.03 ±0.08 ±0.07 ± 

Model 5*; Capture models with fallback (/capture 
16 10 0.03 0.17 1.06 0.76 0. 

±0.00 ±0.01 ±0.03 ±0.04 ± 

Model i?*; Capture models with mixing (/capture - 
16 10 0.03 0.17 1.04 0.76 0. 

±0.00 ±0.00 ±0.01 ±0.01 ± 

Model f : Capture models with mixing (/capture = 
16 10 0.03 0.14 0.96 0.67 0. 

±0.00 ±0.01 ±0.02 ±0.02 ± 

Complete Lateral Mixing 

Model A ^: Capture models with fallback (/ejection ■ 
16 10 0.02 0.11 0.90 0.58 0. 

±0.00 ±0.02 ±0.09 ±0.07 ± 

Model B '. Capture models with fallback (/capture - 
16 10 0.03 0.11 0.86 0.56 0. 

±0.00 ±0.02 ±0.07 ±0.07 ± 

Model E^: Capture models with mixing (/capture = 
16 10 0.03 0.12 0.89 0.60 0. 

±0.00 ±0.01 ±0.05 ±0.04 ± 

Model F^: Capture models with mixing (/capture = 
16 10 0.03 0.12 0.88 0.59 0. 

±0.00 ±0.01 ±0.05 ±0.05 ± 



= 1, /fallback 

97 0.66 
0.07 ±0.06 

= 0, /fallback 

91 0.58 
0.03 ±0.02 

= 1, AMn,ix : 

89 0.48 
0.01 ±0.00 

= 1, AM,,ix = 
81 0.47 
0.02 ±0.01 

= 1, /fallback 

74 0.71 
0.08 ±0.11 

= 0, /fallback 

71 0.63 
0.07 ±0.06 

: 1, AMn,ix = 

74 0.59 
0.04 ±0.04 

: 1, AM^ix = 
73 0.59 
0.05 ±0.04 



= 1, M2 = 1.45 M®, Mbh = 5.4 M®) 
0.43 0.32 0.08 0.08 8.00 2.51 
±0.06 ±0.06 ±0.02 ±0.02 ±2.05 ±0.07 

= 1, M2 = 1.45 Mq, Mbh = 5.4 Mq) 
0.36 0.27 0.08 0.07 5.88 2.44 
±0.01 ±0.01 ±0.01 ±0.01 ±0.85 ±0.03 

= 2 Mfeuback, M2 = 1.45 M® , Mbh = 5.4 M@) 
0.24 0.18 0.06 0.06 6.70 2.43 
±0.00 ±0.00 ±0.00 ±0.00 ±0.13 ±0.02 

= 3 Mf,uback, M2 = 1.45 M@, Mbh = 5.4 M@) 
0.26 0.19 0.06 0.05 10.89 2.45 
±0.01 ±0.01 ±0.00 ±0.00 ±0.65 ±0.03 



= 1, M2 = 1.45 M® 
0.63 0.63 0. 
±0.11 ±0.10 ± 



= 1, M2 = 1.45 M( 
0.55 0.54 
±0.07 ±0.06 

2 Mfaiiback, M2 - 

0.49 0.50 
±0.05 ±0.04 

3 Mfallback, M2 -- 

0.49 0.50 
±0.05 ±0.04 



, Mbh = 5.4 M®) 
.55 0.37 8.25 3.38 
0.07 ±0.07 ±2.03 ±0.64 

Mbh = 5.4 Mq) 
51 0.30 6.64 3.03 
0.05 ±0.04 ±1.34 ±0.48 

.45 M©, Mbh = 5.4 M®) 
51 0.28 7.78 3.43 
0.03 ±0.03 ±1.88 ±0.64 

.45 M®, Mbh = 5.4 M®) 
51 0.28 8.56 3.65 
0.03 ±0.03 ±2.04 ±0.61 



4.83 2.17 
±0.37 ±0.35 



4.36 2.72 
±0.44 ±0.23 



49 90 
±19 



35 66 
±10 



4.40 2.54 3.94 50 53 

±0.16 ±0.04 ±0.32 ±2 

5.03 1.89 7.75 82 90 

±0.13 ±0.07 ±0.41 ±4 



4.91 2.09 

±0.34 ±0.34 

4.62 2.44 

±0.38 ±0.28 



49 90 637665 
±19 

38 90 177557 
±16 



4.69 2.31 2.53 53 90 195131 
±0.37 ±0.32 ±1.52 ±18 



4.78 2.19 3.41 60 90 141353 
±0.35 ±0.33 ±2.33 ±18 



Note. — See Tables 2a and 2b for definitions. 



